Prom Luttinger liquid to Altshuler-Aronov anomaly in multi-channel quantum wires 
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A crossover theory connecting Altshuler-Aronov electron-electron interaction corrections and Lut- 
tinger liquid behavior in quasi-lD disordered conductors has been formulated. Based on an inter- 
acting non-linear sigma model, we compute the tunneling density of states and the interaction 
correction to the conductivity, covering the full crossover. 
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I. INTRODUCTION 



In many regards, the theory of interacting one- 
dimensional (ID) conductors seems to have reached a 
mature state [HQ. In the absence of disorder, inter- 
acting ID conductors are often described as Luttinger 
liquids (LL) [J. The LL model has been shown to be rel- 
evant for a variety of physical systems, e.g. carbon nan- 
otubes, semiconductor quantum wires (QWs), nanowires, 
cold atoms in ID traps, and quantum Hall edge states. 
Thanks to the work of Kane and Fisher Q, a number 
of precursors @, H, @ , and a lot of research activity in 
the past decade, the physical properties of both clean 
LL's and LL's polluted by only few impurities appear 
to be reasonably well understood. In particular, vari- 
ous expected scaling laws have been checked numerically 
using the functional renormalization group Q. In the 
complementary case of significantly disordered systems, 
two distinct scenarios may be realized: in quantum wires 
supporting only a small number N = 0(1) of conducting 
channels, a conspiracy of interactions and impurity scat- 
tering leads to instant localization on microscopic length 
scales ~ I, where I is the (renormalized) elastic mean 
free path. However, in systems with many conducting 
channels — as realized in, for instance, multi-wall car- 
bon nanotubes (MWNTs) — low energy excitations may 
diffuse over distances Nl 3> I much larger than the mean 
free path (before Anderson localization sets in at scales ~ 
Nl.) The low-tempcraturc behavior of transport coeffi- 
cients, tunneling rates, etc., is governed by strong anoma- 
lies differing from both the characteristic power laws of 
the non-Fermi liquid clean LL state, and the strongly 
localized disordered few-channel state. While for thick 
diffusive wires, both the Altshuler-Aronov corrections to 
the conductivity [|[ and the zero-bias anomaly (ZBA) in 
the tunneling density of states (TDoS) [E |T(| are well 
understood, much less is known about the crossover be- 
tween the LL and the diffusion-dominated regime. In 
one exemplary study on the TDoS, Mishchenko, Andreev 
and Glazman (MAG) [ll[ quantitatively worked out the 
crossover from a LL dominated high-temperature regime 
to the diffusive low-energy regime. However, the inter- 
mediate physics of several other important system char- 
acteristics, say, the conductivity, remains unexplored. 



Motivated by the fact that LL/diffusive crossover sce- 
narios are realized in several experimentally relevant sys- 
tems, the purpose of this paper is the construction of 
a comprehensive theory interpolating between the two 
regimes. We will apply this theory to rederive earlier re- 
sults on the TDoS, and to obtain an expression of the 
conductivity that contains LL power laws and Altshuler- 
Aronov corrections as limiting cases. Specifically, we will 
formulate and analyze a model of spinless electrons with 
repulsive forward-scattering Coulomb interactions and a 
short-ranged Gaussian random potential. This model 
will be shown to reproduce known limiting behaviors. 
Of course, in real-life systems, additional mechanisms, 
such as clcctron-phonon scattering, magnetic impurity 
scattering, or electronic (e.g. superconducting, charge- 
density wave, or Peierls) instabilities may also be im- 
portant but we do not address those complications 
here. The theory of dirty /ew-channel Luttinger liquids 
has been pioneered by Giamarchi and Schulz (l3| . see 
also Refs. [H, [IE EE ELI 01 for related work. In con- 
trast, we here consider the case of wide A^-channel wires, 
N 3> 1, termed 'quasi-lD' throughout. We note that 
in the absence of disorder, and despite the presence of 
many channels, LL behavior may still be realized 
This is well-known from studies of many-chain systems 
P, 0, 03 1 where an LL phase is observable as long as 
interchain hopping can be neglected [2(| . 

Besides offering various technical simplifications, the 
multi-channel theory is also of considerable practical in- 
terest. For example, MWNTs typically carry N m 5 
to 20 open channels. Many recent MWNT experiments 
have shown that the interplay of disorder and interactions 
plays a decisive role in these systems [2l|, [U HE HH HE 
|26|. On the theoretical side, the TDoS of a MWNT was 
shown to exhibit a characteristic ZBA at low energies 
fill . I27I l2l| . As far as we know, no results for the inter- 
action correction to the Drude conductivity for arbitrary 
disorder have been published so far. Other examples of 
experimental interest include quasi-lD or gan ic conduc- 
tors [1^1, charge-density wave nanowires |30l|. polymer 
nanofibers [3l|, L12O3 wires [32|, or MoSe nanowires [33[. 
In particular, InAs nanowires currently receive a great 
deal of attention. They have typical mean free paths of 
the order 10 to 100 nm, with N w 10 to several 100 and 
lengths of microns up to millimetres. 
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Both LL power laws and the low-temperature anoma- 
lies in diffusive systems root in the same physical mech- 
anism: scattering off impurities creates &2kp oscillatory 
screening cloud (Friedel oscillation) in the conduction 
electron density which acts as an additional source of 
backscattering. This idea was put forward in Ref. [3J] for 
a weakly interacting ID system with a single barrier. For 
a disordered ituo-dimensional system, the ZBA [35j and 
the conductivity [36|, [37| were explored for arbitrary Ttq 
within a perturbative diagrammatic framework. How- 
ever, in (quasi) ID, the power-law decay exponent of 
the Friedel oscillation itself is affected by interactions, 
and the simple picture outlined above no longer applies. 
Relatedly, Fermi liquid renormalization procedures (3(| 
must break down in the ID setting A coherent de- 
scription of quasi-lD systems rather requires the con- 
struction of a generalized approach, and will be the main 
subject of this paper. We note in passing that recent 
work by Golubev et al. [1^, H(| has also explored similar 
crossover phenomena in granular models. 

Building on a replicated coherent-state field integral, 
we will derive an interacting nonlinear a model (NLctM), 
similar in spirit to that of Ref. [H[ . This field theory will 
allow for a nonperturbative treatment of interaction di- 
agrams that dominate in the ballistic limit, and thereby 
captures the relevant LL physics without invoking stan- 
dard bosonization schemes [4l| . A superimposed pertur- 
bative treatment will then allow us to compute the low- 
temperature diffusive corrections to the conductivity in 
conceptual analogy to earlier diagrammatic approaches. 
(In many respects, the ideas entering this combination 
of perturbative and nonperturbative elements are sim- 
ilar to those developed in the theory of Kamenev and 
Andreev 0.) 

The rest of the paper is organized as follows: In Sec.HIl 
we describe the model and the field-theoretical approach 
taken here. The TDoS is discussed in Sec. IIII1 while the 
linear conductivity will be studied in Sec. IIVI Finally, 
in Sec. [V] wc conclude and offer an outlook. Technical 
details have been delegated to several appendices. Below, 
we put h = Ub = 1- 



II. MODEL AND FIELD-THEORETICAL 
APPROACH 

A. The model 

We consider a disordered and interacting quasi-lD 
quantum wire supporting N conducting channels below 
the Fermi energy. To describe the system we employ the 
prototypical Hamiltonian 



H = Hq + Hd 



(1) 



where 



H = -iJ2 v FC J dx¥d x $, 



represents the kinetic energy, C = ±1 labels right - 
and left-moving states, vf is the Fermi velocity, and 
^ = {^i,c, ■ ■ ■ , ^n,c) is an A^-componcnt vector of field 
operators. Disorder scattering is described by 



Hdis 
Hfs 



Hfs + Hbs, (2) 
t^f>~/~w,_^ ( 3 ) 



= I dx%(x)Vc(x)^ c (x), 
c 



H 



BS 



dx^ R {x)W{x)^ L {x) + h.c. 



(4) 



where Hfs and Hbs represent the forward (FS) and 
the backward (BS) scattering, respectively. Here, Vc — 
{Vnn'(x)} and W = {W nn >(x)} are N x N matrices 
in channel space which we assume to be randomly dis- 
tributed. Specifically, 



~3nn'&mm'fi{x X ), (5) 



{W nm {x)W^ nl {x')) dis = 

ZltV\TO 

where U\ = N/ttvf is the ID density of states, and To 
the (bare) elastic scattering time. The (bare) mean free 
path then is I = vftq. Finally, defining the local electron 

density as p = Y^cPC' PC ~ ^c^c, the interaction 
Hamiltonian is defined as 

Hi = jJ dxp 2 (x), (6) 

where the constant U sets the interaction strength. The 
Hamiltonian implies several idealizing assumptions: lin- 
carizable spectrum, channel-independent Fermi velocity, 
absence of electron-electron backscattering, and a num- 
ber of others more. In Appendix [X] we detail how the 
above Hamiltonian may be derived from a more micro- 
scopic model, e.g. one describing a MWNT. In particu- 
lar, we will argue that at energies e < ej_ smaller than 
the characteristic energy above which structures in the 
direction transverse to the wire are resolved, the above 
model Hamiltonian properly describes all relevant system 
characteristics. 

Following and generalizing Ref. (l3j . a further simpli- 
fication may be effected by a gauge transformation re- 
moving the forward scattering disorder. For N = 1, for- 
ward scattering merely leads to phase multiplying the 
electron wave functions and, therefore, does not affect 
(gauge invariant) observables. A non-abelian generaliza- 
tion of that argument may be used to eliminate the for- 
ward scattering matrices V c . To this end_j we perform the 
unitary transformation ^c(x) = Uc{x)^c(x), where 

X dx'Vc(x')] (7) 



Uc{x) = Vexp 



with the path-ordering operator V defined via 

iC 



■Pexp 



dx'V{x') 



dx r , 



dx r . 



n=0 



dx\ 



xV{x n )V{x n . 1 )---V{x 1 ). 
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This unitary transformation eliminates all FS processes, 
Hq + i/ps — * Hq, when H is expressed in terms of the 
'new' fermions, \L r c- One can check explicitly that un- 
der this unitary transformation, Hbs does not change its 
structure, and W(x) remains a Gaussian random vari- 
able with statistical properties ([5|). Moreover, the gauge 
invariant interaction term Hj, and all observables stud- 
ied below remain unaffected by this transformation [42| . 
From now on, we switch to the fields "J and rename them 
'J. Effectively, we can then forget about the FS processes. 

Before proceeding we note that the clean Hamiltonian 
Hq + Hj describes a (multi-channel) LL phase, described 
by the dimcnsionlcss Luttingcr interaction parameter K 
and the plasmon velocity v, 



K 



1 



v F /K. 



(8) 



This LL phase comes along with N — 1 (effectively non- 
interacting) plasmon excitations moving at velocity vf 
[l9j |. For K = 1, we retrieve a Fermi gas, while the re- 
pulsively interacting case corresponds to K < 1. 



B. Hubbard- Stratonovich transformations 

Our aim is to formulate a tractable low-energy field 
theory from the above model. Restricting ourselves to 
equilibrium phenomena, we will employ an imaginary 
time formalism and the replica trick: To compute the free 
energy F = —T\nZ, we use that \nZ = \im r ->Q(Z r — 1) /r 
and represent the disorder average of the rth power of 
the partition function (Z r )di S in terms of an r-fold repli- 
cated coherent-state Grassmann field ip„Ca{x, T )i where 
a = 1, . . . , r. (As usual, an analytic continuation r — > 
to be performed in the end of the calculation is implied.) 

We begin by decoupling the interaction Hamiltonian 
in a standard way via a Hubbard-Stratonovich (HS) 
transformation. Introducing the (replicated) real field 
ip a (x 1 r), this generates the action for the interaction 
part, 



Si = 



I -^j + iip a p a 



(9) 



Note that the fields ip do not carry a channel index, but 
couple only to the ID density p(x). (For the more general 
case (72 7^ 34 discussed in Appendix [A] one introduces a 
doublet fc and proceeds in a similar way.) 

One advantage of the replicated framework is that the 
disorder average may be conveniently performed at an 
early stage of the calculation. The action describing the 
disorder contribution to (Z r )riit then assumes the form 



S, 



1 



(lis 



E 



dxdrdr' 



X 1p m -C ! oc>(%,T')'lp m -C ! oc(x,T), 



where we used Eq. ([5]). We decouple the time-nonlocal 
action Sdis by another Hubbard-Stratonovich transfor- 
mation to obtain 

Sdis = -^7- /" dxdrdr' ^ Q R , aa > (x, t, t') (10) 

oca' 

x Q L , a ' a {x,r' ,r) + ^- J dxdrdr' ^ 

x tpnCa(x, r)Q c , aa > (x, t, r')^ nCa ' (x, r'). 

This HS transformation introduces a functional integral 
over a pair of time-bilocal but spatially local auxiliary 
matrix fields Qr/l- To warrant convergence of the func- 
tional integral, the condition 



Q\ = Ql 



(11) 



is imposed. Implicitly understood in the above discussion 
is the case of broken time-reversal invariance, e.g., by 
application of a weak magnetic field. Otherwise, proper 
inclusion of the Cooperon channel (describing, e.g., one- 
loop weak localization effects) would require an addi- 
tional doubling of the field degrees of freedom [H[ . 

At this stage, we may integrate out the fermion fields 
to arrive at the action 

S[Q,if] = ±Tr<pU- l <p+^'R{Q R Q L ) (12) 
- JV V Trln (do + W+ -^Qc 

where dc = d T — iCvpd x . The trace symbol indicates 
summation over replica indices and integration over space 
and (imaginary) time. 

In the noninteracting limit U — 0, a standard gradient 
expansion of the tracelog in Eq. (fT2"|) leads to a diffu- 
sive action. This expansion is stabilized by the condition 
N 1 (for N = 1, there is of course no diffusive phase), 
and starts by identifying the saddle-point solution. It 
is convenient to switch from time to energy space via 
a double Fourier transform. The saddle configurations 
are then homogeneous (x-independent), parity invariant 
[Qr = Ql), and time-translational invariant (diagonal 
in energy space). The noninteracting saddle Qrjl = A 
can be expressed using the standard notation [3q 



A aa /(e,e') = sgn(e)(5 ee '(5 aa ' 



(13) 



[Our Fourier convention is ip(x,r) = 

TJ2eI^ eiqX+t£T ' t l J ( ( l' e )' with fcrmionic Matsubara 
frequencies e = (2m + l)nT (integer m).] In order to 
study diffusion properties, we need to include fluctuation 
modes around the saddle A. There is one important 
subtlety at that stage, which we discuss in detail in 
Appendix [Bj namely there are two different types of 
massive fluctuation modes. These are (i) excitations 
with Q 2 1, and (ii) chirally asymmetric excitations 
with Q 2 = 1 but Ql 7^ Q R [involving the field W\ in 
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Sec. Ill Ej . Excitations of type (i) are not coupled to 
massless fluctuations at the Gaussian level, and hence 
can safely be neglected. However, type-(ii) excitations 
are linearly coupled to massless excitations [the field 
Wq in Sec. Iff E| . and thus cannot simply be ignored. 
In fact, only when integrating them out, we do obtain 
the physically meaningful ID result for the longitudinal 
diffusion constant, D = v^tq. Moreover, we also recover 
the transversal diffusion constant. Including only 
Gaussian fluctuations around the saddle (flU|) . we thus 
reproduce the full anisotropic (cf. remark in the end of 
Appendix [Bj diffusive theory [ll|, [2?| , with the correct 
longitudinal and transversal diffusion constants. This 
provides an a posteriori justification of the simplifying 
assumptions underlying the effective Hamiltonian (fl]). 



C. Chiral anomaly 

We now turn to the interacting case and look for 
the saddle-point configurations corresponding to the ac- 
tion (Ti~2|) . Unfortunately, we are unable to determine the 
exact saddle, except for the clean limit (tq — > oo). Hence 
we will adopt an Ansatz, motivated by Ref. [9j], which 
will lead to results generally exact to lowest order in U, 
and to any order in the ballistic limit. (In the diffu- 
sive case, the approximate saddle point scheme is good 
enough to address the strongest low-temperature singu- 
larities in transport coefficients.) Specifically, we will, for 
fixed ip a , consider field configurations of the form 

Q c (x,t,t') =e lKc(x - r) A(r-r') e ~ lKc{xy \ (14) 

i.e. configurations differing from the nonintcracting sad- 
dle Qc = A only by a chiral gauge transformation, me- 
diated by a space-time local field Kc = {Kca{x, r)}. 
Varying the action with respect to Kc determines the 
particular Kc that solves the saddle-point equation to 
lowest order in the interaction. 

Before proceeding to the saddle-point solution, let us 
analyze the effect of a gauge transformation breaking 
chiral symmetry. To remove the phase factors dressing 
A in Eq. (TT4")) . we apply a gauge transformation to the 
fermionic fields corresponding to the tracelog in Eq. (fl2"|) . 



steps, and then adding the corresponding Jacobian con- 
tributions [il] . One finds that the Qc fields, describing 
low-energy properties, do not contribute to the anomaly 
(which is a high-energy feature). For the sake of simplic- 
ity, however, we here show a different derivation solely 
relying on known results for the clean case. 

We first recapitulate the 'standard' form of the 
anomaly [43| . 

-NTrln(d c +i<p) = -JVTrln(0 c ) + ^(^c^p), (17) 

where the frequency-momentum representation of the 
chiral polarization bubble is given by 



2 vpq + iCwn 



(18) 



and uj m = 2itmT (integer m) are bosonic Matsubara fre- 
quencies. For later convenience, we also define 



U(q,u m ) = ^nc(?,w ro ) = Vi 



{vFqf 



{v F q) 2 + uj'I 



(19) 



We then use the identity 



Trln ( dc + vp+—Q c ) = 



Trln(3 c + i<p) + Trln ( 1 + -^r(dc + t» _1 Qc 



Trln(<9 c + itp) + Trln ( 1 + -^r(dc + ir c ) _1 A 



(20) 



where Qc is given by (JT4J), we used 



-iKr 



(dc + itp)-^** = (dc + iTc)' 1 , 



and the cyclic invariance of the trace. (Unlike with the 
UV singular logarithmic contribution to the trace, us- 
age of 'cyclic invariance' is a permissible operation in the 
expansion terms.) Similarly, 



Tr In ( d c + iT c + — A 

4t 



(15) 



This transformation has two effects, namely (i) the term 
(p in the tracelog is replaced by Tc, where 



c 



(fi- 



de K c , 



(16) 



and (ii) the classical chiral U(l) invariance is broken at 
the quantum level, leading to a chiral anomaly [43l. [44|. 
The anomaly arises from the Jacobian of the measure 
of the fermionic path integral under the chiral gauge 
transformation (JTSJ) , and is at the basis of path-integral 
bosonization. The calculation can be done by separat- 
ing gauge transformations like Eq. ()15|) into infinitesimal 



Trln(9 c + *To) + Trln ( 1 + ^r(<9c + iT c y L A 



(21) 



Subtracting Eqs. J20]) and IpTjl. and using Eq. (JTTJ) , we 
obtain 



iVTr In ( d c + if + e iKc Ae~ lKc 

4t 



NTrln [d c + iT c + -r- A 
4r n 



^Tr((fU C f) - iTr(r c n c r c ) 



(22) 
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The last two terms constitute the chiral anomaly con- 
tribution to the action. Summing over C, the anomaly 
reads S a = ^Tr(ipH a ip), with 



n a = n-V^n c ^. 



(23) 



Clearly, both Tc and Il a depend on the field Kc- As 
a result, Eq. (fT4")l together with the gauge transforma- 
tion ((TSJ) puts the action (fl"2"l) into the equivalent form 



S 



1 



Tr[tp(U- 1 +Il a )<p] 



(24) 



^Tr [e^Ae-^A] 
8tq 



-N 



c 



Trln 9, 



tT 



c 



4r n 



-A 



where K? = K 



R 



Kr 



D. Saddle-point solution 

Our so far discussion applied to arbitrary time- 
dependent gauge transformations Kc- We now wish to 
find values of Kc such that the transformed Q-ficld 1(14"]) 
represents an (approximate) saddle point of the theory. 
To this end, wc define 



Qc e 



-iK c 



(25) 



and determine Kc such that at Qc = A the action is ap- 
proximately stationary. Put differently, we require that 
at the saddle-point values of Kc, first-order variations 
Qc — > A + 5Qc vanish. Substituting A — > A + 8Qc into 
Eq. (|24p and expanding to linear order in 5Qc, we thus 
obtain the saddle-point equations (one for each value of 
C=±) 



D -iCK 2 (x,T) 



A(7 



t '\^CK 2 (x,t') 



(26) 



2iN f i t 

<9 C + iTc + -. — A 

7Wi V 4t 



These equations can generally not be solved exactly. Fol- 
lowing Ref. [H wc will thus be content with determining 
a solution that holds to leading order in the interaction 
strength U . Expanding Eq. (|26|) to lowest order in if and 
hence U, we find 

[K 2 ,A]- = (Gc(v + d c K c )Gc) XTT ,. (27) 

Here we define the noninteracting chiral disordered Green 
function as [43 



G^ 1 = d c + (V4r )A. 
One may then check that 



(28) 



K Ca {q,uJ m ) 



iiv m + qv F C + i sgn(^ m )/T 
w ™ + (?«f) 2 + |w m |/r 



(29) 



solves Eq. ([27)) . In the clean (ballistic) limit, the solution 
reduces to 



Kc(q, ui m ) 



CvFq 



ILJ„ 



Indeed, these functions represent exact (to all orders in 
U) solutions of the saddle point equation (|2"6"|) . In the 
complementary diffusive limit |w m |To, |<?|/ <C 1, K R / L 
asymptote to @] 

«sgn(w m ) 
L><7 2 + |cj m | 

where D = vf^ro and the left-right asymmetry vanishes. 
With Eq. (29]), T c follows from Eq. $T6j) as 

, . -MC(cJ + sgn(w)/2r ) 

r c (g,w) = 2^ F g — J— — vw,w), 



w 2 + (v F q) 2 + \u\/t 



and Eq. (J23J yields 



-3(^g) 2 + (M + l/r f 
(w 2 + ( UFg ) 2 + M/t ) 2 



(30) 



(31) 



E. Interacting non-linear sigma model 



Let us now substitute Eq. 
tion (fl2l) to obtain 



1 



into the original ac- 



S 



-Tr [<p(U- l + U a )<p] 



(32) 



7TI/1 , 

8^ 



Tr [e iK *Q R e- iK *Q L ] 



-N V Trln(% + ir c +— Q c 



To make further progress, we need to include fluctuations 
in Qc around the saddle A of the interacting system. 
These fluctuations may be conveniently parameterized 
as 

Qc = TcATj 1 , T c = cxp(-Wo/2) (*p(iCWi/2), 

(33) 

where [Wo,i,A]+ = for the relevant low-energy modes. 
(Fluctuations departing from the manifold (|33[) are mas- 
sive and decoupled from massless excitations at the Gaus- 
sian level. Therefore they have no significant effect as 
long as N > 1.) The condition [W ,i,A] + = is re- 
solved by setting 



Wo. 





-Bl 



Bo,i 




(34) 



where the 2x2 structure refers to energy space (e > 
0,e < 0). The soft fluctuations of the theory are thus 
parametrized by unconstrained complex-valued matrix 
elements (B 0:1 ) ee >. aa >(x), with e > 0, e' < 0. Note 
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that this parametrization gives Q* R = Ql as required 
by Eq. (HH), and Q\ = Q\ = 1. 

At this point, it is convenient to integrate out high- 
energy fluctuations of the field (p within the framework 
of a standard renormalization group (RG) analysis [l3[. 
Starting with the bandwidth cu c (which is of the order of 
the Fermi energy) as an UV cutoff, we subsequently inte- 
grate out all ip fields in the strip lo' c < \uj m \ < lu c . Here, 
the intermediate cutoff uj' c (tq -1 ,! 1 ), but small enough 
to allow for a perturbative treatment of the remaining 
(low-energy) interaction field. We will show in Sec. IfV Al 
that the precise value of tu' c drops out in the calculation 
of observables like the interaction correction to the con- 
ductivity. To implement the RG, we temporarily switch 
to yet another gauge representation of the theory, viz. 

S = lTr[( /3 (f/- 1 +n)^] + ^Tr( e ^Q fl e-^g L )(35) 
2 «r 



N Tr In 



'C 



4t 



e i(Kc-K° ) Qce -i(K -K° ) 



differing from Eq. (132p by a gauge transformation medi- 
ated by exp[i(K c — K c )), where dcK c = -(p. Within 
the high-energy strip, ^-fluctuations, coupling to the ac- 
tion through a 'low-energy' perturbation of 0(1/ To), are 
expected to be of no importance, and will be disregarded 
here. It then turns out that high-energy tp fluctuations 
renormalize both To's, the one in the second term and 
that inside the tracelog of Eq. (|35|) . in the same manner. 
Since Kc — K c = K-i for high energies, this renormaliza- 
tion follows by evaluation of (e^ 2 ^' 1 "' -7 ''" 2 ^' )]) for en- 
ergies within the strip, using the first term of Eq. (|35p to 
carry out this Gaussian average. As net result, we end up 
with an effective renormalization of the mean free time. 



r 



r {uj' c /u} c y% 



(36) 



with the Luttingcr liquid power-law exponent for weak 
backscattering by a single impurity @ , 



7, = 2(1 - K)/N, 



(37) 



where K is the Luttinger parameter ([8]). In what follows, 
we imagine that the RG has been carried out, and only 
discuss its effects explicitly where required. (In Sec. IIII1 
it will in fact be more convenient not to perform the RG 
transformation at all.) 

Its logarithmic nonlinearity makes the general analysis 
of the action (|3"2"|) difficult. For momenta q < 1/(vfTq), 



and frequencies ui < r _ , probing the multiple scattering 
regime, a gradient expansion in small momenta can be 
performed to obtain 

S = ilr [ip (JJ- 1 + U a )ip] (38) 



2 

7TZ/1 



^Tr \T r AT b 
8t L r 



~ 1 e~ iK ' 2 T L AT L ~ 1 e lI<2 



Y Tr (A(T^ 1 d c T c + iT c )) 



c 



+ ^Y. T ^ G ^ T c 1 dcTc + l Tc)f + .- 



where the Green functions are given by (|28[) . In Ap- 
pendix [B] we show how the standard form of the diffu- 
sive (T-model may be obtained from this prototype action. 
However, for arbitrary momenta and frequencies no such 
reduction is possible. Nonetheless progress can be made, 
viz. by noting that an expansion of the action (|32|) to 
second order in the fluctuation generators B is sufficient 
to compute our observables of prime interest. In essence, 
this expansion is tantamount to an RPA approximation 
to the model where the 'RPA-bubble' accounts for impu- 
rity scattering. Our neglect of higher-order fluctuations 
essentially implies that the renormalization of diffusion 
and interaction by quantum interference processes is not 
taken into account. 



III. ZERO-BIAS ANOMALY 

To warm up, before addressing the more involved in- 
teraction corrections to the conductivity in Sec. IIV1 let 
us briefly discuss the energy-dependent bulk tunneling 
density of states (TDoS), see also Refs. [TJ H3, The 
TDoS governs temperature- or voltage-dependent trans- 
port quantities when contacts or other conductors are 
weakly coupled to the QW. It is well known that in 
both limits, Luttinger and diffusive, the DoS exhibits pro- 
nounced structure as a function of energy. The crossover 
between the two limiting scenarios was addressed by 
MAG [IH by a direct analysis of the screened interac- 
tion in multi-channel quantum wires. We here rederive 
and extend their results, essentially to illustrate the ap- 
plication of the present formalism. 

For a particle-hole symmetric system, the TDoS is [l[ 

/>oo 

p{oJ,T) = -(2/7r)coth(w/2T) V / dtsmut ImC/ C (*), 

c ^° 

(39) 

where Gc(t) is the real time analytic continuation, r — > 
?i + 0+, of 

^(T T il> nC (x,T)$ nC (x,0)) = (40) 

n 

T T the time-ordering operator, and x an (arbitrary) space 
coordinate. Since Qr = Ql, we are free to focus on Qr. In 
Eq. (|40|) . fluctuations Tc ^ I do not contribute, and the 
average in Eq. (j40|) is taken using the action S v obtained 
from Eq. (J32]) with Tq = 1. (The replica limit is then 
trivial to take.) For noninteracting fermions, Kc = 0, 
and hence 

Gait) = —iTrviA(t) = —iftisiT/sinh(TrTt), 

, , s f°° sinM) 
p{u,T) = 2nTcoth(w/2T) / dt . ) ' =v x . 

J smh(7rTf) 

In order to analytically evaluate ( e l l K R( x > T )- K R( x >°)]} 
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we expand S v to quadratic order in ip, leading to 
Stn = -Tr ipU cp. 



(41) 



This truncation does not introduce any approximation 
in the ballistic limit, and is consistent with the weak- 
interaction condition inherent in the construction of the 
saddle fT4"|), The kernel n includes (i) the term (U^ 1 + 
H a ), (ii) a contribution from the second term in Eq. (f3"2")h 
and (iii) from the expansion of the tracelog, a term ~ 
Tr(GcTc) 2 - For instance, (ii) gives 



TTVl , 



^Tr(AK 2 AK 2 - K 2 ) 

OT 



= Tr ip(-q,-u)- 



v 1 {v F q) 2 \Lo\/2TQ 



(tu 2 + {v F q) 2 + M/t ) 2 
Combining all three terms, the kernel in Eq. (|4"Tj) reads 



110(9,^ 



U 



{v F qf 



+ (vFq) 2 + |w m |/r 



(42) 



This expression bridges the ballistic and diffusive 
regimes. It is exact in the LL limit, while it describes 
diffusively screened interactions in the opposite limit, as- 
suming weak interactions. 

It is then straightforward to arrive at the TDoS. With 
the definition (f8|) of the LL parameter K, and the bulk 
1 7^ and boundary (7b) Luttingcr tunneling exponents 

j=(K+l/K-2)/2N, lb = {l/K-l)/N : (43) 
some algebra leads to [46[ , 
p{<o,T) 



2 f dt 

-coth(w/2T)Rc / — smiujt) 

* Jo t 

( rrTt 



(44) 



Vsinh(7rT£) 



7+1 



(1 + iuj c t)~ 



7_Fi(t/T ,T)+7 t ffi(t/-ro,T)+i7F 2 (t/T )+i 7!) ff 2 (i/-ro) 



x e 



where we have defined 

wt T\ T3 r, ( l 1 \ l-cos(2y M ) 

Fliy > T) = RC Jo ^ [u - T^TTu) tanh( M /(2Tr o) ) 



Hi(y,T) 
H2(y) 



2 
iin 



du- 



1 — cos(2yu) 



u 2 + iu y/2 tanh(u/(2Tr )) 



-ye-y (h(y) + I (y)) . 



and Io.i arc Bcssel functions of imaginary argument [47lj . 
For T = 0, the integrals for F\ and Hi can be evaluated 
in closed form, 



Fid/,0) = -In 



2e 



-c 



+ K (y) cosh(y), 



#i(y,0) = l-y[cosh(y)A^)+sinh(y)^ (y)], 



where C = 0.577 ... is the Euler constant, and Kq 1 
are related Bessel functions of imaginary argument [471 ] . 
Note again that Eq. (|44|) has been derived under the as- 
sumption of large channel number TV. 

Before proceeding, let us pause to briefly comment 
on the relation of Eq. flU} to the results of MAG. The 
far reaching equivalence between Eq. (|4"4")l and MAG's 
Eq. (14) is best exposed by shifting their integration vari- 
able t according to t — > t — i/2T. The results then turn 
out to be identical, provided one identifies the integral 
kernel V(w) of MAG with 



V(w) 



1/K- 



l Re w(l - K) + i/t 



2N uj 3 / 2 ^uj + i/t 



Effectively, however, in MAG the limiting case of strong 
interactions and/or large channel number, K — > 0, is con- 
sidered. In the ballistic regime, and for moderate inter- 
actions and channel numbers, their results do, therefore, 
not scale according to the Luttinger exponent Eq. (|43[) . 
Within the framework of MAG's analysis, this deviation 
is however straightforward to avoid by not taking the 
limit of strong interactions at an early stage. 

We next discuss Eq. (|4"4"|) in various limits of interest. 
Beginning with the zero temperature case, T = 0, for 
ujtq S> 1, only small t contribute to the integral (|44|) . 
The functions F i:2 and if^a are then basically negligible, 
the i-integral can be done, and we obtain the well-known 
LL power-law TDoS Q, 



r(i + 7) 



(45) 



with the Gamma function T(x). In the diffusive regime, 
ujt <C 1, large t dominate in Eq. (|4"4")) . and one can use 
asymptotic expansions for F± t 2 and H\ 2 - Actually, the 
precise crossover scale separating the diffusive from the 
ballistic regime in the T = TDoS is set by uitq ~ 7, see 
Ref. flu. With the error function $, we find 



1 - $ 



2tt 



LOTq 



76 



e C 7 + 7 i, 

(2w c t )t 

(2W C T ) 7 7T7b 



(46) 



corresponding to an exponentially vanishing TDoS at low 
energies. This pseudogap behavior was first note d by 
Nazarov [l(j , and subsequently rederived in Refs. [111. [281]. 
While a direct expansion of the s'mlujt) in Eq. (|4"4"|) would 
suggest a linearly vanishing TDoS [281 ] , the correct result 
is the nonperturbative exponential law (|46[) . In fact, one 
can show by complex integration methods that in the 
ID case the prefactor of any polynomial-in-w term must 
vanish. Similarly, one can show that the T-dependence 
of p(uj = 0) also has to be exponential as T — » 0. Figure 
[1] shows the full T = crossover solution together with 
the asymptotic results (|4"5"j) and (TlGf) . 

Next we turn to the regime Ttq <C 1 and ujtq <C 1. The 
crossover function describing the transition from u> < T 
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FIG. 1: Zero-temperature crossover solution (|44p for the 
TDoS (full curve) at K = 0.3 and N = 5. Here lj c t = 5 x 10 4 
was chosen, which is large enough to render results insensi- 
tive to its value. Dotted and dashed curves give the limit- 
ing expressions (|46|l and (|45|l . respectively. Note the double- 
logarithmic scales. 



to id > T, reported in Ref. [Ill, can be reproduced from 
Eq. (HH) using asympotic expansions of and iJi,2- In 
particular, for u) <C T -C 1/tq, 



C7+7b 



( Tt 



(2w c t ) 7 V 2.389 \2ir-yl 



1/4 



-1. 075 7i , y/2ir/Tr 

(47) 

Finally, once cither Tto ^> 1 or lutq 3> 1, typical values 
for t in Eq. ((44]) are t ~ min(l/T, 1/w), such that effec- 
tively the standard LL power laws are recovered. For 
T» (I/td.w), 



p(T) /TTX7 r(l + 7 ) 



r 2 (i + 7 /2) 



(T/w c 



(48) 



IV. CONDUCTIVITY 
A. Interaction correction 

In this section, we discuss the interaction correction to 
the Drudc conductivity for arbitrary disorder, connecting 
the diffusive Altshuler-Aronov result Q with the ballistic 
limit, where Luttinger liquid power laws emerge. It is 
worth reiterating that weak localization corrections are 
not included in our theory, see Sec. Ill Bl 

The Kubo formula determines the linear dc conductiv- 
ity as the q — > 0, iuj m — > u> + i0 + limit of 

€ 2 C 2 

a(q, u) m ) = (j(q,uj m )j{-q, -u m )) H , (49) 

us m mu> m 

where po is the Drude resistivity, m the effective mass, 
and j the particle current. The conductivity can alter- 
natively be expressed as 



a(q,uj m ) = -(e/q) 2 uj m P(q,uj m ), 



(50) 



where the density-density correlator can be written as 



((p(q,uj)ip(-q,-uj)) _ 1_ 
U 2 U' 



(51) 



Starting from the full action in the original formulation, 
see Eq. (|12[) . the formula (f5"Tj) may be derived as follows: 
Introduce a source field ry(g, u) m ) that is added to tp in the 
tracclog only. The density correlator then follows from 
the resulting generating functional Z[rj\ by means of 



d 2 \nZ 



dr}(q, w m )dr)(-q, -u m ) J 



(52) 



ri=0 



Under the iocai gauge transformation (|15[) . rj remains 
invariant, and one can safely perform the shift <p ^ p — ?l 
in the action. After that shift, it is a simple matter to 
carry out the derivatives in Eq. (|52[) and to arrive at 
Eq. (f5Tj) . In what follows, we use the intermediate cutoff 
u>' c and the mean free time t' q = T( ) {uj' c /uj c ) li with the 
exponent (|37p . see Sec. Ill El After the corresponding RG 
procedure, the remaining interactions are treated by a 
first-order perturbation theory approach. 

The classical Drudc conductivity can now be recovered 
by our formalism as follows. Dropping all B fluctuations 
in Eq. (|3"8"|). the action up to quadratic order in tp is S v , 
given by Eqs. (jlTj) and (|12")) . Keeping only S v , we may 
then apply Eq. (|50[) . continue to real frequencies and ar- 
rive at the Drude conductivity 



Ne 2 VF 



(53) 



This result includes an interaction correction encoded in 
the high-frequency RG scheme above. When substitut- 
ing Tq — > To, one recovers the standard (noninteracting) 
Drude conductivity, 



e 2 v 1 v 2 F T . 



We mention in passing that the above approximation also 
obtains the correct ac Drude form of the optical conduc- 
tivity 

<td{v) = vd/0- — iwro). (54) 

Interaction corrections to the dc Drude conductivity 
(|53p . defined via a = a' D + 8a' , are obtained by inclusion 
of (i) B fluctuations, and (ii) higher-order contributions 
in tp to the action. To lowest order in the interactions, 
these additional ingredients to the theory conspire to ob- 
tain the result 



-7*/(l+7<) 



l + 7 4 ln(2 7 re- 1 ^ c T^) 



7;Rc / dn 



n 



n 



Q + i/r' Q 



2T^(i + K) v / n 2 + in/T[ ) 



(55) 



with the function 



9(ficoth(Q/2T)) 

on ■ 



(56) 



For weak interactions, it is justified to use the renormal- 
ized mean free time in the form 



^ = T (a, c To)-^ +1 ). 



(57) 



Before turning to the derivation of this result, let us dis- 
cuss its physical meaning. We first note that, up to the 
prefactor (w c To) _7i ^ 7i+1 ^, the conductivity is a univer- 
sal function of the interaction parameter K, the chan- 
nel number N, and the parameter Ttq. Importantly, 
the crossover temperature separating the LL and the 
Altshulcr-Aronov regions, respectively, is determined by 
the renormalized scale Tq. For Ttq > 1 and 7 < 1 
Eq. (|55|) simplifies to 



o(T) 



(w c t£)-* [l + l % ln(2^ e -( 1+c >T^) 



T>t'-\ 
(58) 

Exponentiation of the logarithm obtains the familiar Lut- 
tinger liquid high-temperature power law 0, Il3l . [l8j , 



o(T) 

<TD 



a -7i(l+C) 



{2-kT/u 



T>t' 



(59) 



otherwise obtained by stopping the RG procedure of 
Sec. Ill El at T » t' -1 . Governed by the single- impurity 
backscattering exponent 7$ (cf. Eq. ([37)) and Ref. Q), the 
exponentiated form is valid only at asymptotically large 
temperatures where multiple interference is negligible. In 
the complementary regime Ttq <C 1 the interaction cor- 
rection becomes 



CD 



-li 



-1/2) 



1 + K 



-1/2 



(60) 



where £ is the Riemann zeta function [47| . This ex- 
pression recovers the well-known T -1 / 2 Altshuler-Aronov 
correction in ID Q , leading to a pronounced suppression 
of the conductivity at low temperatures. At very low 
temperatures, this correction becomes sizeable and our 
first-order perturbative approach ceases to be valid. In 
that case, also quantum interference corrections have to 
be included. 

The crossover solution ([55)) smoothly interpolates be- 
tween the limits (|59[) and ([60)) . respectively, as is shown 
for K = 0.9 and N = 30 in FigureEl Note that Eqs. ([59]) 
and ([6"0")) match each other at Tt' q of order unity. Re- 
lated crossover scenarios from diffusive dynamics at low 
temperatures to non-diffusive dynamics at high temper- 
ature have been discussed in earlier work. Specifically, 
in Refs. [3^, the conductance of a quantum dot ar- 
ray was considered. In that system the high-temperature 
regime is governed by the chaotic dynamics of individ- 
ual quantum dots. The latter differs strikingly from the 
ballistic dynamics of clean quantum wires which is why 




for 

m 



FIG. 2: Temperature dependence of the conductivity ([55 
K = 0.9 and N = 30. The limiting behaviors ([59)) and 
are shown as dotted and dashed curves, respectively. 



the results obtained for the high-temperature conduc- 
tivity corrections are quantitatively different from ours. 
(Nonetheless, the qualitative high-T profile looks simi- 
lar, cf. Fig. 2 of Ref. [3|.) 

We finally note that our calculation above ignores 
finite-size effects, which are of minor importance in suf- 
ficiently long wires with good contact to external leads, 
where the resistance is dominated by the intrinsic part. 
Note also that four-terminal resistance measurements al- 
low to circumvent boundary effects as long as the elec- 
trodes act as non-invasive probes. In both cases, the 
conductivity is the relevant transport coefficient. This 
regime has been experimentally studied in carbon nan- 
otubes and other quantum wires for many years by now. 
In particular, linear voltage drops along a MWNT can 
be measured and have been reported in Ref. [50l | . 



B. Derivation 

Let us now sketch the derivation of Eq. ([55)) . For nota- 
tional simplicity, the primes in Tq and uj' c will be omitted 
for the moment. Our strategy is as follows: We expand 
Eq. ([32]) to second order in ^-fluctuations, and integrate 
them out. This will generate an effective action for ip, 
which we expand up to fourth order. In addition, for in- 
ternal consistency, the ^-independent terms in Eq. ([3"2"]) 
must also be kept up to fourth order. Schematically, our 
starting action then reads 



S = Sip + Sb + S, 



(i) 



(61) 



The interaction correction results from the last three 
terms. The contribution Sb combines (i) terms of second 
order in Bq i but zeroth order in tp, and (ii) of first order 
in B but second order in <p. Note that the first order in 
both B and <p is strictly zero, due to our choice of the 
saddle point. Moreover, contributions to second order in 
B and first or second order in ip identically vanish in the 
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replica limit. Finally. 



,(1) _ TWxVF 



^Tr(G c r c ) 4 



(62) 



comes from S-independent terms in the expansion of the 
tracelog in Eq. (|32|). while 



Sf ] = ^ (^TrAl - l -TrKK 2 KKl + iTrAA'f Aiff 

(63) 

similarly follows from the second term in Eq. (|32l) . 

Since Eq. (|5ip involves only ip, the ^-fluctuations can 
be integrated out at the level of the action, replacing 
Sb effectively by a new fourth-order term S\ [<p]. The 
lengthy but straightforward derivation of this term is de- 



- c(l) 



j(2) 



,Ci) 



tailed in Appendix Id Defining S4 = S 
the lowest-order correction to the conductivity. Eq. (fSTj) . 
is obtained by first-order expansion in S4, 

SP(q, uj) = -U- 2 (<p(q,w)<p(-q, -uj)[S 4 - (S*)s v ]) Sv . 

(64) 

Loosely speaking, fluctuations of the tp field (S^\s^) 
describe the RPA screening of the Coulomb interaction 
by impurity-dressed Green functions ('empty bubble'), 
while joint fluctuations of the B and the ip field (S4 ) ac- 
count for the vertex renormalization of the RPA bubbles 
by impurity ladders. Accordingly, we will find that S4 
holds responsible for the corrections to the conductiv- 
ity showing the strongest low-temperature singularities. 
Conversely, the terms dominantly contribute to 

the conductivity in the Luttinger limit Ttq > 1. 

Doing the ^-contractions, we obtain for the conductiv- 
ity 



5a 



a D lm J ^F(fi) J glLj-^Q) (65) 
i (K 1 + K a +K 3 ){q,Sl) 



(vFq) 



(v 2 F q 2 - " 2 - ifiAfc) 



with the auxiliary quantities 



Ki(q,Q) = 4(1 
K 2 (q,Ct) = -2, 
K 3 (q, fi) 



1 - 2iflT 



(2v F qT Q ) 2 + (1 - 2ifi 7fc) 2 ) 
1 



T 2 (v 2 F q 2 -n 2 -in/T ) 
2 - 6i^r n - (2fLr ) 2 + (2v F qT ) 2 
(1 - 2iQr ) 2 + (2v F qT ) 2 



(66) 



where the function F(£l) has been defined in (|56"]) . The 



contributions A"i 2 3 stem from S, 



(1,2,3) 



respectively, and 



Il R denotes the retarded function corresponding to 



IIq , see Eq. (|42|) . Specifically, the 'ballistic contribu- 
tions' Ki 2 result from summing 12 terms corresponding 
to the diagrams in Fig. [31 followed by analytic continua- 
tion to real frequencies. (Although standard in principle, 




4x 



o{+) 



4x 



FIG. 3: Diagrams corresponding to and in the con- 
ductivity calculation. The contraction of S^* 1 or with 
<p{q, u>)ip(— q, — u>) (denoted by </?(±), respectively) leads to 
the pairing of <p(±) with two different Tc (or K2 in the 
case of S^ 2 '), implying 12 different contractions. Due to the 
cyclic invariance of the trace, only three different contribu- 
tions emerge, depending on the Tc ordering. Solid curves 
represent schematically the trace, external dotted lines the 
two Tc connected to </?(±), and dashed-dotted lines the pair- 
ing between the remaining Le- 



the actual calculation obtaining these terms, in particu- 
lar the analytic continuation to real frequencies for the 
product of four Green functions, is rather involved. We 
refer to Appendix A of Ref. [36| for related technical de- 
tails.) The calculation of the third term A3 is detailed 
in Appendix [Cl Dominating in the diffusive limit, it ob- 
tains the standard ID Altshuler-Aronov correction 
see Eq. ([BO]). 

Equation ((65]) implies the total correction to the Drude 
conductivity to lowest order in the interaction, 



5a 



-2a D U Im 



F(n) 



dq 
2^ 



(67) 



v 2 F q 2 [v 2 F q 2 + (1/tq - lilf 



(v%q 2 -n 2 - iVL/T ) 2 {v 2 q 2 -ft 2 - ifl/ro)' 

where v — v F /K, see Eq. (|8]). After carrying out the final 
momentum integration and restoring primed quantities, 
we arrive at the preliminary result 



Sa' 



-7icr^Re 



cm 



(68) 



2T' {i + K)^jn 2 + in/T' Q 



which may give the misleading impression that the con- 
ductivity depends on the somewhat arbitrary cutoff u>' c . 
Noting that lo' c t' q > 1, to first order in the interaction, we 
have 

T = 7 i 1 + H T Wc)) : 

leading to the estimate ([57)1 . Splitting off the UV diver- 
gent part in Eq. (|68[) leads to the logarithmic term in 
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Eq. (|55[) , while the upper limit can be sent to infinity for 
the remaining integral in Eq. (|68p . We finally arrive at 
the ^-independent form (|55p . 

V. CONCLUSION 

In this paper, we have developed a low-energy field the- 
ory of weakly interacting disordered multi-channel con- 
ductors. The theory is formulated in terms of two auxil- 
iary fields, a scalar field tp(x, r) decoupling the electron- 
electron interactions, and a matrix field Qc{x,t,t') de- 
coupling the effective interactions arising from the dis- 
order ensemble average. The two fields tp and Qc are 
coupled by a gauge mechanism. For general values of 
interaction/disorder/channel number, the theory is gov- 
erned by a nonlinear action — the notorious 'tracelog' — 
and remains difficult to evaluate. However, in a number 
of important cases analytical progress is possible. Specif- 
ically, at energies e < t^ 1 smaller than the elastic scat- 
tering rate, and for large channel numbers, we recover 
the diffusive interacting a model [38( with all the known 
consequences. In the opposite limit e 3> t^ 1 , an (asymp- 
totically exact) mapping onto the familiar action of the 
Luttinger liquid is possible. Finally, for weak interac- 
tions, a low-order expansion of the tracelog in ip and in 
the generators of Q-fluctuations becomes permissible. 

Focusing on the latter regime, we apply the formalism 
to study the crossover of the tunneling DoS (previously 
described in Ref. [ll[) an d the conductivity from the bal- 
listic to the diffusive regime. In the ballistic limit, we 
recover the standard Luttinger liquid power law behav- 
iors in these quantities, while in the diffusive limit, we 
obtain a pseudogap in the low-energy part of the TDoS 
and the T -1 / 2 Altshuler-Aronov interaction correction. 
Our main results are Eq. (|44[) for the tunneling density 
of states, and Eq. for the temperature dependence 
of the conductivity. Wc believe that these results, cov- 
ering the entire crossover from ballistic to diffusive, will 
be valuable in understanding experimental data on nan- 
otubes or nanowires. 

Let us conclude by noting some of the open prob- 
lems in this area. The question of what happens to the 
conductivity at very low T cannot be answered by our 
lowest-order calculation. Relatedly, in order to treat the 
small- A limit, new conceptual advances will be necessary. 
(Progress along this line has been reported in Ref.[l8[.) 
Other interesting open questions not addressed here in- 
clude the weak localization correction and the magneto- 
conductivity, the inclusion of the spin degree of freedom, 
and a microscopic calculation of the dephasing time in 
ID. These topics may be addressed by future work. 
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APPENDIX A: DERIVATION OF THE 
EFFECTIVE MODEL HAMILTONIAN 

In this appendix we review how the model Hamilto- 
nian ([I]) may be distilled from more microscopic descrip- 
tions of a quasi-lD conductor. Generally, the simplify- 
ing assumptions below are expected to hold for energies 
6 e i- 1 where e± is the characteristic energy scale related 
to transversal (to the cross section of the wire) excitations 
in the system. For I < L±, we have e± ~ vf/L±, where 
L± is a typical transverse width of the wire. For I > L±, 
on the other hand, e± ~ vpljL\. Resolving phenom- 
ena on larger energy scales is a doable task which, how- 
ever, requires a more refined (and less universal) model- 
ing. Throughout, the phrase 'low energy' refers to the 
regime e < e±. 

In the absence of disorder and interactions, a quasi-lD 
conductor is characterized by N open channels energeti- 
cally below the Fermi energy. The number N sa S±/A F 
is often tunable, e.g., via doping or backgate voltages, 
where S± = L\. We assume that (i) N 3> 1 is suffi- 
ciently large to justify the approximations employed later 
on, and (ii) the Fermi energy is located sufficiently far 
away from the bottom of the bands formed by longitu- 
dinal momentum components to justify introduction of 
well-defined chiral (right- or left-moving, C = R/L = ±) 
electron branches for each channel; clearly, this becomes 
problematic if the Fermi energy is close to the bottom of 
a band. Moreover, we shall (iii) focus on the case of spin- 
less (or spin-polarized) electrons. The generalization to 
the spinful case does not pose conceptual problems and 
is left to future work. 

Denoting the ID coordinate as x (where we assume 
that no confinement along this axis is present) and the 
transverse degrees of freedom as rj_ , the electron operator 
^ e (r), with r = (x,r±), can be expanded into the iV" 
transverse eigenfunctions as 

N 

*e( r HE E e iCknX M*x)iPnc(x), (Al) 

n=l C=R/L 

where 4> n (r±) are the transverse eigenmodes normalized 
according to 

J dr ±( j)* n {r ± )(j) m (r ± ) = S nm . (A2) 

To give an example, for MWNTs, the eigenmodes, aris- 
ing from the wrapping of a 2D graphene sheet onto a 
cylinder, are [27[ 

— 1/2 

4> n (y)=L ± ' exp(2Triny/L±), (A3) 

where L± = 2itRq for outermost-shell radius Rq of the 
MWNT. Here, the transverse coordinate is the angular 
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variable y, with < y < L±. For given confinement, 
each band intersects the Fermi surface at k = ik n with 
its own Fermi momentum k n , where the slope is given by 
the respective Fermi velocity v n . At low energy scales, 
the kinetic part of the Hamiltonian then is 

Ho = -iJ2 C [ dx^ c (x)vd x ^ c (x), (A4) 
c •* 

with v — diag(ux, v 2 , ■ ■ ■ , vn). The noninteracting DoS is 
thus given by V\ = J2 n =i l/( 7r ' y n)- To simplify matters, 
we neglect differences in the channel- dependent velocity, 
Vn — * (vn)n = vf, where vf denotes the average channel 
velocity. As may be checked, e.g., by an explicit calcu- 
lation of the diffusive two-point correlation function of 
the system, this assumption is permissible at low ener- 
gies. For equal Fermi velocities, the kinetic part of the 
Hamiltonian assumes the form of Hq in Eq. {]]) . 

Next let us turn to the repulsive Coulomb interac- 
tions among the electrons in the QW. Starting from an 
arbitrary microscopic interaction potential ?7 m j cr (r, r'), 
which incorporates external screening effects due to the 
substrate or surrounding gates, one arrives at a rather 
complicated ID interaction Hamiltonian. At low ener- 
gies, it is however sufficient to consider a simple model 
interaction, which is assumed to be long-ranged (e.g., a 
1/r potential) on length scales larger than L±. In that 
case, the long-range tail of the interaction is expected 
to dominate all relevant ID Coulomb interaction matrix 
elements, such as 

U nm (x - x') = J dr ± dr' ± U mlcr (x - x', r ± , r' ± ) 

x |0„(r ± )| 2 |0 m (rl)| 2 . 

For \x — x'\ 3> Lx, the potential U m i cr is basically in- 
dependent of the transverse coordinates rj_ and r^_, im- 
plying that the projected interaction does not depend on 
the channel indices, U mn — > U(x — x'). The expansion 
(|A1[) then leads to an effective contribution to the Hamil- 
tonian, 

N 

ft = - ^ / dxdx'U(x-x') (A5) 

n,m=l Ci,C 2 ,C 3 ,C 4 
x g— t(Ci — d)k n x ^-i{Ci— C'3)k m x' 

X V'lcx ( x )^Lc 3 ( x ')^mC 3 (x')lpnCt (x). 

Following the standard reasoning leading to LL models 
P, 0, for a long-ranged potential, electron-electron (e- 
e) backscattering (C\ = — C4, C2 = — C3) processes will 
be strongly suppressed. Moreover, e-e Umklapp scatter- 
ing is important only close to commcnsurabilitics and 
will not be discussed here. For a long-ranged interaction, 
the dominant interaction process then corresponds to e- 
c forward scattering, where C\ = C4 and C2 = C3 in 
Eq. (|A5p . Such processes describe interactions coupling 
the ID density fluctuations. Hi = 5 J dxdx' p(x)U (x — 



x')p(x'). It is then justified to express Hi in terms of the 
q = Fourier component of U(x — x'), denoted by U . (In 
the clean case, up to multiplicative logarithmic correc- 
tions in most observables, this procedure also describes 
the case of the unscreened 1/r interaction [l|.) We are 
then left with the simple ID forward-scattering interac- 
tion ([|5]). We mention in passing that e-e backscattering 
can be (partially) included by allowing for different in- 
teraction strengths 174 and (72 in the pcPc an d pcP-c 
couplings, respectively [l[. With minor modifications, 
our approach can be adjusted to this situation. 

Finally, quenched disorder is included by starting from 
a short-ranged (3D) Gaussian random potential V(r), 
with the disorder average defined by its only nonvanish- 
ing cumulant 

(V(r)V(r')) dls = -^—S(r - r'), (A6) 

Z1TV 3 Tq 

where v% = v\/S±_ is the 3D noninteracting DoS and we 
assume that disorder is weak, k n l 3> 1 for all k n . With 
Eq. (|A1[) . we then obtain FS and BS scattering as in 
Eqs. ([3]) and respectively, where 

Vc,nw = e- iC ( fe »- fe »')x / ' dvMr±)<t>n'(ri-)V(r), 

W nn , = e -*Cfc-+*nO» f drMr ± )ct> n ,(r ± )V(r) 

can be taken as independent random variables whose sta- 
tistical properties follow from Eq. (|A6[) . To simplify mat- 
ters, we assume the specific torus- type confinement (|A3|) . 
Unless one is interested in details of transverse fluctua- 
tions or other high-energy features, the precise form of 
the confinement potential is not expected to change the 
essential physics. Statistical correlations of W nm (x) can 
then be described by the simple form ([5]). 

APPENDIX B: NONINTERACTING LIMIT OF 
THE NLctM 

Let us here discuss the diffusive properties at \q\l <C 
1, |w m |To < 1 of the NLcrM ([38)1 in the noninteracting 
limit. Using the paramctrization (|33p . we have two dif- 
ferent types of quantum fluctuations, Bq and B±, where 
the latter break chiral symmetry. The B\ modes have a 
gap, but nevertheless are not to be discarded since they 
couple to the massless Bq fluctuations. We show this now 
on the Gaussian level, starting from the noninteracting 
version of Eq. ([55]). 

T ° c 
+2CA(T c 1 d x T c ) - 2eQ c ] , (Bl) 

where we put Vp — 1 in intermediate steps. In order to 
recover the conventional NLcrM, we define the chirally 



symmetric field Q = TqAT^ 1 with T = exp(-W /2), 
and integrate over the massive modes B\. Expanding S 
to second order in B\ (but zeroth order in Bo), the first 
term in Eq. (|B1[) produces 



with kernels 



see Eq. (|34|) , which determines the mass gap of the fluc- 
tuation mode B\ referred to in Sec. Ill Bl The remaining 
parts of the action (|B1[) are then expanded to linear order 
in B\ , where only the Tq 1 8 x Tc term gives a contribution. 
Here one has to be careful because of the non-commuting 
nature of B _i, but within a first-order expansion, no dif- 
ficulties arise under the parametrization (|33|) . Terms that 
are quadratic in B\ or involve spatial derivatives of B\ 
are neglected since they vanish at low energy and/or long 
length scales. We finally arrive at 



Tr(D(d x Q ) 2 - 2eQ - W 2 /2t 



TTV\ 

Wi[T - 1 a,T ,A] 



IIi(g,w) 

n 2 (g,w) 
n 3 (g,w) 



Moreover, 



9 2 



+ \w\ {\u>\ + l/2r ) 



v 2 F q 2 



(M + 1/2t ) 2 ' 
qv F sgn(u})/T a 

(|a;| + l/ro)(|^| + l/2To)+4g 2 
v 2 F q* + (M + l/2r ) 2 



Next Bi is integrated out, and we finally obtain the ID 
diffusive action for Qo, 



S = ^j- (DTr(d x Q ) 2 - 4Tr(eQ )) . (B2) 



On the Gaussian level, keeping replica indices implicit, 
this results in 



x r c (<7o - 9,^0 - u> m )ic(q,u m ,qo,uo) 

+ K 2 {q 1 uJ m )K 2 {qo -q,u> - cJ m )sgn(e„ - u m ) 



dq 
2^ 



e„>0,e m <0 ' 

x [Dq 2 + (e„ - e m )] |S 0; „ m ((7)| 2 , 



and 



implying the usual diffusion pole. 

We finally remark that at energies e > ej_, the dif- 
fusion properties (e.g. of MWNTs) become anisotropic. 
These 'medium energy effects' may be resolved by em- 
ploying the full disorder correlator of the W nm (x), and 
decoupling the impurity four-fermion action in terms of 
a channel-dependent Q-field. One thus obtains a more 
structured theory wherein details of anisotropic trans- 
port are resolved. 



Pi(q ,e n ,e n - v ) = ^ (-iC) / — \T c (q, 

C,u> m ^ 

x r c (<7o -q,wo- u m )Ic{q,u m ,qo,uo)), 



APPENDIX C: INTERACTION CORRECTION 

(3) 

Here we provide the derivation of S\ and of the cor- 
responding interaction correction, see Sec. IIV Bl The ex- 
pansion of Eq. ([52")) in B yields 



with the kernel 



S B = ^TMWiIWi + WqIWo I c (q,u; m ,q ,LjQ)=iv F [ ^G c (p,e) 

8t J 2tt 

-W Ii 2 W 1 + 2W P + 2W 1 P 1 ], xG c (p - go, e - u )G c (p -q,e- u m ) 
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The Gaussian ^-fluctuations are then integrated out. 

(3) 

Calling the resulting action S± , we find 



o(3) _ ™i 
64 ~ 4^ 



E E E 

90,91,92 wo>0,ui,iiJ2 0<e„<w o 

y^rc(gi,o;i)rc(go - gi, lu - wi) 

c 

iC n 2 (g ^o) 



X-fc(9o,w ,gi,Wi) 1 - 

V 4 n 3 ((?o,wo) 

+^2(91,^1)^2(90 - 9i, w - wi)sgn(e„ - wi) 



^r C /(<7 2 , w 2 )rc7'(-9o - 92, -w - w 2 ) 
L c> 

iC> n 2 (g ,w ) 



Xlc> (9o, ^o,9o + 92, w + w 2 ) 1- . jj t 
+-K 2 (92,w 2 )if 2 (-g - 92,-^0 - w 2 )sgn(e„ - u - lu 2 ) 
x (lii(g ,cj ) - n^(go,wo)/4n 3 ( g o,wo))~ 1 



^(-jC , )r c (gi, ^i)r c (9o - 91, w - Wi) 
L C 



x7c(9o,w ,9i,wi) 
xr C '(9o - 92, w - w 2 ) 
x/ C ' (9o, w , 90 + 92,^0 + w 2 ) 



^(-iC')rc(92,w 2 ) 



/n 3 (g , w ) 



Equation (|64|) is then used to compute the conductiv- 
ity correction. The contraction of ip(q,Lu)ip(—q, — u>) with 



(3) 

S\ ' leads to eight terms, with four different terms ap- 
pearing twice each. These four terms are in fact two by 
two equivalent via the symmetry lj — > —lu. In particular, 
the two contributions with —lo vanish after the summa- 
tion over e„ is performed. The two possible contractions 
are (i) 91 = —92 = 9, Wi = — lo-i = w, and (ii) 91 = 9, 
lui = lo, 92 = 9 — 90, w 2 = lo — loq. The external momen- 
tum q is now taken to zero, and we also let lo — > under 
the constraint Rew > 0. The e n -summation can be done 
easily and yields, after the frequency shift luq — > wo + w m , 



8vj?e 2 lSlT§ ^ LJ + LJ m - \lO - LJ m \ 

6a 3 = ^ " 



U)„ 



n (9o, w ) 



-^1(90, w ) (111 (g , w ) - Il^go, Wo)/4n 3 (go, w )) 1 

+ ^2(9o,^o)/n 3 (9o, w ) 

where Lo m — ► has been anticipated and we have used 
the abbreviations 



1 



Xi(q ,(j ) = — g- 



Crc(go,w ) 
4tq ^ 2iv F qo<p{qo,uo) 



( iCn 2 (go,w )\ o « X 2 (9o,w ) 

x 1 ~ TTTi \ p c 9o,w + — r. 

V 4 n 3 (g ,w )/ 2t iUFgo<A9o,^o) 

v / \ 1 (-^ r c(go>wo) r o / \ 
A 2 (g , loo) = -7-0 ) 77— — -7— — r7 c (9o,w ), 



4t o 2 4r' 2wf9o</>(9o,w ) 



Ic (9o, w ) = 2r 



1/tq + w - iCv F q 
(l/2r + w - iCv F qoT z 



The summation over cjo is then replaced by a contour 
integral, with a line cut for £1 = iloq = ilo. Performing 
the analytic continuation ilo = to' + i0 + with a/ — > 0, we 
finally obtain K3 in Eq. (p5|) . 
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